SUMS OF MAGNETIC EIGENVALUES ARE MAXIMAL 
ON ROTATIONALLY SYMMETRIC DOMAINS 



RICHARD S. LAUGESEN, JIAN LIANG AND ARINDAM ROY 

Abstract. The sum of the first n > 1 energy levels of the planar 
Laplacian with constant magnetic field of given total flux is shown to be 
maximal among triangles for the equilateral triangle, under normaliza- 
tion of the ratio (moment of inertia) /(area)'' on the domain. The result 
holds for both Dirichlet and Neumann boundary conditions, with an 
analogue for Robin (or de Gennes) boundary conditions too. 

The square similarly maximizes the eigenvalue sum among parallelo- 
grams, and the disk maximizes among ellipses. More generally, a domain 
with rotational symmetry will maximize the magnetic eigenvalue sum 
among all linear images of that domain. 

These results are new even for the ground state energy (n = 1). 



1. Introduction 

Eigenvalues of the Laplacian on a plane domain represent energy levels 
of a quantum particle in two dimensions. These eigenvalues are known in 
closed form only for special types of domain, such as disks and rectangles. 
Consequently a great deal of effort has gone into proving upper and lower 
bounds on eigenvalues in terms of geometric properties of the domain, such 
as area and perimeter. 

For example, the Rayleigh-Faber-Krahn inequality says that the ground 
state energy Ai of a quantum particle with Planck constant h confined to a 
region of area A is bounded below according to 

where jo,i is the first zero of the Bessel function Jq. Equality holds for 
the disk. Such bounds provide not only hard estimates, but also soft in- 
tuition, for they indicate how geometric attributes of the domain constrain 
the analytic information encoded in the eigenvalues. 

A magnetic field imposed transversely through the domain makes the 
energy levels even more difficult to determine theoretically. In this paper we 
aim to discover geometrically sharp estimates on such magnetic eigenvalues. 

Given a bounded plane domain O, write Xj ($7, h, /?) for the jth eigenvalue 
of the Dirichlet Laplacian with Planck constant h > 0, under a constant 
transverse magnetic field (0,0, /3). These Dirichlet eigenvalues are defined 
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Figure 1 . A plane domain with 3- fold rotational symmetry, 
and its image under a linear map T. 

precisely in the next section, as are the Neumann eigenvalues (assuming 
the domain has Lipschitz boundary). The eigenvalue equations are 

Dirichlet: < ^ ' (1.1) 

I u = on dVt 

Neumann: <^ ^ , . (1-2 

[n ■ {hV - iF)u = on on ^ ^ 

where the vector potential F{x) = f (— X2, xi) creates a field V x F = (0, 0, /3). 
Let A be the area and / the moment of inertia of about its centroid: 

\x — cp dx 

with centroid c = f^xdx. 

Our method works by linearly transforming a bounded plane domain D 
having A^-fold rotational symmetry, such as in Figure [TJ The main result. 
Theorem 13. H says that sums of magnetic eigenvalues decrease under linear 
transformation of the rotationally symmetric domain, provided the Planck 
constant and field strength are transformed suitably too. This theorem 
encompasses both Dirichlet and Neumann boundary conditions, and Robin- 
de Gennes conditions are treated in Theorem 13. 2[ The ability to handle all 
three boundary conditions is a striking feature of our method. 

By expressing our main theorem geometrically, we will obtain: 

Corollary 1.1. Assume D is a bounded plane domain with rotational sym- 
metry of order greater than or equal to 3, and that Q is the image of D 
under a linear transformation. Fix /? E M. 

Then for each n > 1, the normalized eigenvalue sum 

x,{Q,h,^) + --- + Xn{n,h,^)]^ (1.3) 

is maximal when Q. = D, for each n>l. 

Maximality holds also for dilated, rotated, and reflected images of D; and 
when n = 1, every maximizing Q is of that type. 

The same conclusions hold for sums of Neumann eigenvalues. 
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The corollary is proved in Section [71 

Notice the total flux of the magnetic field is the same for each domain, in 
the corollary, since multiplying the field strength (3/ A by area A gives (3. 

In particular, the corollary implies that the normalized eigenvalue sum 
()1.3p is maximal among triangles for the equilateral, maximal among paral- 
lelograms for the square, and maximal among ellipses for the disk. To apply 
this result in practice, one would like explicit formulas for the eigenvalues 
of the equilateral triangle and the square — but we do not know whether 
such formulas exist for /3 7^ 0. 

Incidentally, for triangles the moment of inertia can be calculated in terms 
of the side lengths as / = (Zj + Z| + /|)^/36, while for a parallelogram with 
adjacent side lengths li, I2, the moment of inertia equals I = {if + l2)A/12. 

Two reasons for studying eigenvalue sums are that the sum represents the 
energy needed to fill the lowest n states under the Pauli exclusion principle, 
and that summability methods improve the behavior of high eigenvalues, 
which are difficult to study directly. An example of the improvement pro- 
vided by summation is that although the Polya conjecture remains open, 
claiming the Weyl asymptotic provides a lower bound for each eigenvalue of 
the Dirichlet Laplacian, Li and Yau [25] were able to show that sums of these 
eigenvalues are indeed bounded below by the analogous Weyl asymptotic. 

Prior work. The topic originates with an old result of Polya on the fun- 
damental tone of a membrane (as stated in [27] and proved in |28[ Chapter 
IV]). He established the case n = 1 of this paper, for vanishing magnetic 
field; that is, he studied the first eigenvalue of the usual Dirichlet Laplacian. 
His method relied heavily on uniqueness of the fundamental mode, and thus 
it failed to extend to Neumann boundary conditions, higher eigenvalue sums 
or the magnetic spectrum, all of which we treat in this paper. 

Our work relies on the method of Rotations and Tight Frames, developed 
recently by Laugesen and Siudeja [23] for finding upper bounds on eigenvalue 
sums of the Laplacian. We extend this technique to handle magnetic fields. 
The resulting theorems are new even when n = 1. 

The magnetic situation presents two new challenges. First, the interac- 
tions of the momentum (or gradient) operator with the magnetic field must 
be averaged over the iV-fold rotations of D; see the quantities Q2 and Q3 in 
Section [H Second, a new proof must be found for the equality statement (see 
Section [6]) , because the eigenfunction u is complex- valued in the magnetic 
situation, and so u and u are non-equal in formulas (]6.4p and (j6.5p . 

The only prior geometrically sharp estimate we know, for magnetic eigen- 
values, is the Faber-Krahn type lower bound of L. Erdos That result 
says the normalized ground state energy Ai(r2, h, ^)A is minimal for the disk 
among all two dimensional domains. Notice the total fiux is once again the 
same, for each domain considered. 

The lower bound of Erdos relates to our upper bound in Corollary 11.11 as 
follows. One can rewrite our functional in (II. 3p as (Ai + • • • + Xn)A ■ A"^ /I. 
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Thus it arises from multiplying a Faber-Krahn term (Ai + • • • + A„)j4 that is 
normahzed by area with a purely geometric, scale-invariant term jl that 
penalizes long, thin domains. Erdos studied solely the Faber-Krahn term, 
and only for n = 1, although his work does apply to arbitrary domains. We 
would like to extend our work to arbitrary domains too, but the task seems 
challenging even in the absence of magnetic field [23^ §4]. 

The scarcity of isoperimetric type inequalities for magnetic eigenvalues 
contrasts with a profusion of results in the nonmagnetic setting, where one 
may consult the surveys by Ashbaugh and Benguria [2j or Benguria and 
Linde [4j, and the monographs of Bandle |3J, Henrot |15] . Kawohl [16], Ke- 
savan [1?] and Polya-Szego [29j. In particular, triangular domains have been 
much studied [B [HI HH ED EH |23l [Ml EI] . Our aim in this paper is to begin 
developing a magnetic isoperimetric theory of comparable richness. 

Asymptotically sharp inequalities (semi-classical constants). Our 

work in this paper is "geometrically sharp" , since an extremal domain exists 
for each fixed n. The Li-Yau inequality mentioned above (for eigenvalue 
sums in the absence of magnetic field) is not geometrically sharp, but is 
asymptotically sharp since equality holds for each domain in the limit as 
n — 7- oo, by the Weyl asymptotic. Other asymptotically sharp bounds on 
eigenvalue sum functionals were studied recently by Dolbeault, Geisinger, 
Harrell, Hermi, Kroger, Laptev, Loss, and Weidl (see ^ [I2l [HI [18] and 
references therein). 

Fewer such estimates are known with magnetic field, and some natural 
conjectures turn out to be false. Notably, the magnetic Polya conjecture has 
been disproved by Frank, Loss and Weidl who constructed counterex- 
amples even for square domains. A magnetic Li-Yau inequality nonetheless 
holds true for eigenvalue sums, by work of Erdos, Loss and Vougalter [8] 
extending results of Laptev and Weidl [l9] . 

Higher dimensions. We hope to extend our inequalities to higher dimen- 
sions in a later paper. The results will be more complicated, for two reasons. 
First, the moment of inertia must be evaluated on an "inverse" domain, as 
seen already with vanishing magnetic field p5]. Second, the field is no longer 
perpendicular to the domain. 

No sharp Faber-Krahn lower bound is known for the ground state energy 
of the magnetic Laplacian, in higher dimensions. The lower bound of Erdos 
in the plane remains to be extended. 

2. Assumptions and notation 

Vector potentials. Throughout the paper we consider a constant magnetic 
field on the plane, written 

i? = (0,0,/3) 
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for some fixed /3 G M. The field can be generated from the curl of a suitable 
vector potential F £ C^iR^^R"^), as 

B = V X F = (0,0,91^2 - d2Fi). 

Most commonly we employ the potential 

Pi^) = ^(-2^2, 2:1), 

but other choices are allowed too, such as /3(— X2,0) or (3(0, xi). This non- 
uniqueness of the potential illustrates the principle of gauge invariance, 
whereby the magnetic field is unchanged by adding a gradient vector to 
the potential, because the curl of a gradient equals zero. 

Eigenvalues. Consider a bounded plane domain il.. Denote by Xj{il.,h,/3) 
the Dirichlet eigenvalues of the Laplacian on il. with Planck constant h > 
and constant magnetic field B = (0,0, /?). These eigenvalues form an 
increasing sequence 

< Ai < A2 < A3 < ;> 00 

and are determined from the Rayleigh quotient 

L\(ihV + F)u\^ dx 

R\u\ = , , ,^ ; ' where u G HUn; C , -u ^ 0, 

i^\u\^dx 

by the usual minimax variational principles. In particular, Ai = m.m.uR[u\- 
Note the gradient operator V and the vector potential F are regarded as 
row vectors. 

The existence of these eigenvalues, and of a corresponding orthonormal 
basis of smooth eigenfunctions, follows from standard elliptic theory [301 
Corollary III. 7. 8]. The key facts are that imbeds compactly into LP' and 
the Rayleigh quotient can be written R[u] = Q[u,u]/{u,u)i2(^q-^ where the 
sesquilinear form 

Q[n, v] = {{ihV + F)u, {ihV + i^)w>i2(n) 

satisfies the following three properties. It is: (i) conjugate symmetric, (ii) 
continuous as a function of u, i" G Hq{U;C), and (iii) ff^-coercive after 
adding a multiple of the L^-inner product (meaning + ci||u||^2 ^ 

C2||n||^i for all u G H^{Q,;C), for some ci,C2 > 0). 

The eigenfunction equation (II. ip follows from the Euler-Lagrange condi- 
tion for a critical point of the Rayleigh quotient. Note the eigenfunctions 
depend on the choice of vector potential, but the eigenvalues do not: they 
depend merely on the field strength parameter /3, by gauge invariance (see 
Lemma lA.2p . Thus our notation A(il,/i, /?) need only indicate the depen- 
dence of the eigenvalue on the field parameter /3, not the vector potential. 

The first eigenvalue is positive by diamagnetic comparison with the first 
eigenvalue of the Dirichlet Laplacian (/3 = 0), or else by Lemma lA. 81 (Also 
Xi{Q,,h,l3) > Ai(M^,/i, /3) = h\f3\ by domain monotonicity, but we will not 
need this fact.) 
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The Neumann eigenvalues h, (3) arise from the same Rayleigh quo- 

tient R[u\, but the class of trial functions is larger, namely u G H^{Q; C); we 
make the standing assumption that the domain has Lipschitz boundary so 
that imbeds compactly into L^. The first Neumann eigenvalue /xi is pos- 
itive whenever /3 7^ 0, by Lemma [A. 8 1 Note the boundary conditions in (jl.2p 
arise naturally in the Neumann case, from the Euler-Lagrange condition. 

The magnetic eigenvalues are invariant with respect to rigid motions and 
dilations of the domain. See Appendix lAl for precise statements. 

Matrix notation. Given a real or complex matrix M, write its Hilbert- 
Schmidt norm as 

\\M\\HS = {Y.\M,k\'f' = itrMM^f/' 

where denotes the complex conjugate of the transpose matrix. 

For example, the identity matrix has Hilbert-Schmidt norm equal to \/2- 

3. The main results: sharp upper bounds on eigenvalue sums 

Fix the Planck constant h > and field parameter /3 S M, throughout 
this section. 

Dirichlet and Neumann eigenvalues. We start by linearly transforming 
a rotationally symmetric plane domain D. The Planck constant and field 
strength are transformed suitably too, in order to obtain a sharp result. 

Theorem 3.1. If the bounded plane domain D has rotational symmetry of 
order greater than or equal to 3, then 

y X ( T{D), h, —^-4^) < y A,- (D, h, p) 

for each n > 1 and each invertible linear transformation T o/M^. 

Equality holds if T is a scalar multiple of an orthogonal transformation. 
For nonzero magnetic fields (f3 7^ 0), equality holds for the first eigenvalue 
(n = 1) if and only if T is a scalar multiple of an orthogonal transformation. 

The same results hold for Neumann eigenvalues. 

The proof is in Sections O and El 

The rotationally symmetric domain D in the theorem need not be convex, 
or a regular polygon, or have any axis of symmetry. For example, it could 
be shaped like a three-bladed propeller. 

The equality statement for the first eigenvalue under zero magnetic field 
is more complicated than in the theorem, because rectangles (which do not 
possess rotational symmetry) can also be extremal |23i Theorem 3.1]. On 
the other hand, the equality case in this paper requires new ideas too, in 
order to handle the magnetic interaction terms. 
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Robin-de Gennes: 



Robin— de Gennes eigenvalues. Analogous results can be established for 
the boundary condition of the third kind. This boundary condition was stud- 
ied in thermodynamics by Robin and in superconductivity by de Gennes. 
First we need some definitions. Denote the Robin-de Gennes eigenvalues 
by pj{Q,h, I3,a) where the constant o" > is the Robin parameter. The 
Rayleigh quotient is 

r , fcliihV + F)u\'^ dx + a Lf^lul'^ ds ^, 
R\u\ = ^ ' ,^ , ^-^^ where u G i/^ 17; C , 0. 

(3.1) 

As in the Neumann case, we assume the domain has Lipschitz boundary so 
that the spectrum is well defined. 

The eigenvalue equation and boundary conditions are easily deduced from 
the variational characterization of the eigenvalues. They are: 

{ {ihV + Ff'u = pu in il, 
1 n • — iF)u + (Tu = on 50. 

Notice the Robin-de Gennes case reduces to Neumann when a = 0, and 

reduces formally to the Dirichlet case when o" = oo. 

The first eigenvalue pi is positive for each value of /3, by Lemma lA. 81 
The Robin-de Gennes eigenvalue sums are bounded sharply by the next 

theorem. 

Theorem 3.2. If D has rotational symmetry of order greater than or equal 
to 3, then 

y p ( T{D), h, „ a) < y p, (D, h, (3, a) 

for each n > 1 and each invertible linear transformation T ofM?. 

Equality holds if T is a scalar multiple of an orthogonal transformation. 
For nonzero magnetic fields (j5 ^ Oj, equality holds for the first eigenvalue 
(n = 1) if and only if T is a scalar multiple of an orthogonal transformation. 

The proof is in Section [H) 

One can also express the theorem in geometric terms, similarly to Corol- 
lary [TTTJ (See the nonmagnetic version in j23 i Corollary 3.4].) 

4. Tight frame identities 

When proving the main theorem, we will need to average certain matrices 
with respect to conjugation by the rotation group of order A^. Let the matrix 
Um represent rotation of the plane by angle 2TTm/N. Write Id for the identity 
matrix. 

Lemma 4.1. If N >3, then for every 2x2 real matrix M one has 

1 ^ 1 1 

- UmMUl = (- tr M)ld + -(M - Mt). 

m=l 
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Proof. Begin by decomposing M into its symmetric and antisymmetric parts 
as M = Ms + Ma, where Ms = {M + Mt)/2 and Ma = (M - Mt)/2. To 
handle the symmetric part we argue as for Schur's Lemma in representation 
theory: let 

1 ^ 

m=l 

and observe that V is symmetric. Notice UiVUl = V for each m (because 
UiUm = Um+i and Un+i = Ui). Thus V commutes with the rotation Ui. 
Let a G M be an eigenvalue of the symmetric matrix V, with eigenvector x. 
Then 

V{Uix) = Ui{Vx) = a{Uix). 

Hence Uix is also an eigenvector belonging to a. Since x and Uix are 
linearly independent (noting that Ui rotates by angle ^ < vr, because N > 
3), by taking linear combinations we deduce that every vector in is an 
eigenvector of V with eigenvalue a. Hence V = aid. Taking the trace yields 
a = I tr y. Also 

1 ^ 

tr y = — ^ tr = tr = tr M. 

m=l 

Therefore V = i(trM)ld, which gives the first part of the formula in the 
Lemma. 

Since M is real, its antisymmetric part can be written Ma = ^ ( ? ) for 
some 6 S R. Thus Ma consists of rotation by 7r/2 followed by rescaling by 
h. Hence Ma commutes with the rotation and so 

1 ^ 1 

- UmMaUl = Ma = -{M- M^), 

m=l 

which proves the second term of the formula in the Lemma. □ 

Next we collect some immediate consequences of Lemma l4. 11 

Lemma 4.2. Let N > 3, and suppose M and T are 2x2 real matrices, 
with T invertible. Then 

^f:Urr.T-'T-^Ui = l\\T-rHsl'i, 
m=l 

1 ^ 1 

- U^T^M^MTUl = -||Mr|||,5ld. 

m=l 

// tr M = then 

1 ^ 1 

- Y UmT^^MTUi = -(T~^MT-{T~^MT)^). 

m=l 
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Also, by applying Lemma l4.ll with M = T'^T and then putting on the 
left and y on the right, we find: 

Lemma 4.3. Let N > 3 and suppose T is a 2 x 2 real matrix and y ^ M? 
is a unit column vector. Then 

m=l 



Connection to equiangular tight frames. Fix a vector y £M? of length 
1. Then M = yy'^ is a symmetric 2x2 matrix. Lemma l4. II implies that 



N 



Conjugating with x G yields 



1 ^ 1 

— E x^Umyy^Ul^x = 

m=l 

or in other words 

^ N 
E |x-[/™y|2 = -|x|2, xGR2. 

m=l 

This Plancherel type identity says that the system of vectors {Umy}m=i 
(consisting of the A^th roots of unity, rotated to start at y) forms a tight 
frame with constant N/2. 

More information on frames and their applications in Hilbert spaces may 
be found in the monograph by Christensen [5j or the text by Han et al. [13]. 

5. Proof of Theorem I3.lt the inequality 

We prove the Dirichlet inequality in the theorem. The Neumann proof is 
virtually identical. 

Without loss of generality, we may assume T is diagonal with 



T 



ti 

t2 



where ti, t2 > 0, by a singular value decomposition of T and using invariance 
of the spectrum under rotations and reflections (Appendix |A|) . 
Choose the vector potential to be 

F(x) = ^(-X2,xi) = /3(Mx)t where ^ = ^ 

so that the magnetic field is V x F = (0, 0, /3) as required. (Aside. This 
vector potential points tangentially to circles centered at the origin, and so 
it interacts well with rotations of the domain.) 
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For brevity in what follows, we write 



V2 2 V2 1 



\\T-^Hs ^tf + q^' ||r-i||H5|detr| Y^i + ^i' 

Then be = 2/ [tit^^ + t^^t2) . 

We will need the Rayleigh-Poincare variational principle [3, p. 98], which 
characterizes the sum of the first n > 1 eigenvalues as 

n 

n 

= min I RV^j] '■ vi, . . . ,Vn G Hq(Q; C) are pairwise L^-orthogonal|. 

Specifically, we will construct trial functions Vj on the domain T{D) by 
linearly transplanting eigenfunctions Uj of D, and then we will average with 
respect to the rotations of D by means of the tight frame lemmas in Section^ 
Let ui, ti2, ti3, • • • be orthonormal eigenfunctions on D for the vector po- 
tential F, corresponding to the eigenvalues Xj{D, h, 13), j = 1, 2, 3, . . .. Con- 
sider an orthogonal 2x2 matrix U that fixes D, so that U{D) = D. Define 
trial functions 

Vj = Uj o U o 

on the domain E = T{D), noting Vj G Hq{E;C) because Uj G Hq{D;C). 
The functions Vj are pairwise orthogonal, since 

/ VjVkdx = / UjUkdx ■ I detTC/^^l = 
J E Jd 

when j 7^ k. Thus by the Rayleigh-Poincare principle, we have 

" . . L\{ibhV + cF)vA'^dx , ^ 

,=1 5=1 iE\"]\ "-^ 



For each function v = Vj, we change variable in the right side of (15. ip with 
X I—)- TU~^x and use that the Uj are normalized in L^, thus finding 

fj^\{ibhV + cF)v\'^dx 
|f p dx 

_ \ibhVu{x)UT^^ + cl3{MTU-^xyu{x)f dx ■ \ detTU-^\ 
~ f^\u\^ dx-\ det TU-^ 

= [ \ibhVu{x)UT-^ + cl3{MTU'^x)^u{x)f dx 
Jd 

{Qi + Q2 + Qs) dx 

D 
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by expanding the square, where 

= b'^h'^Vu[UT-^T-^U^{Vu)\ 

Q2 = 2bchf3Re{iuVu[UT-^MTU^x}, 

Qs = c2/32|u|2xt[c/rtM'^MrC/'^]x. 

Let > 3 be the order of rotational symmetry of D, and choose U to 
be the matrix Um representing rotation by angle 2irm/N. Averaging over 
m = 1, . . . ,N implies by Lemma 14.21 that 

N 



m=l 



1 ^ 

_ ^ g2 = bchpR(i{iuVu[T-^MT - {T-^MT)'^]x] 



N 

m=l 



N 



2hf3 Re{m(Vn)Mx}, 



iV ^ 2 

m=l 

a2\ |2| |2 

where we used for Q2 that tr M = and where we simplified the three 
expressions by substituting the definitions of the matrices M and T and the 
constants b and c. 

Hence by averaging (|5.ip over m = 1, . . . , we deduce that 

n 



2| |2 



dx 



" /■ r 1 
< V / /i^lVujf + 2h/3Re {iuj{Vuj)Mx} + -/^Vil 

= V/ \{ihV + ^{Mxy)uj\'^ dx 

n 

which completes the proof of the inequality in the theorem. 

6. Proof of Theorem I3.lt the case of equahty 

We continue to treat the Dirichlet eigenvalues. The Neumann case pro- 
ceeds exactly the same way. 
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Sufficient conditions for equality. When T is an orthogonal matrix, equahty 
holds in the theorem because the eigenvalues are invariant under orthogo- 
nal transformations (see Appendix [Aj) and = has Hilbert-Schmidt 
norm 

(trrtT)i/2 = (trId)i/2 = ^2. Scalar multiples preserve equality too, 
because 

Xj{rD,rh,/3/r)=XjiD,h,/3) 

by Lemma lA.71 

Necessary conditions for equality, when n = 1. Now assume /3 7^ and that 
equality holds in the theorem for the first eigenvalue, so that 

Ai (t{D), h, / ) = X^ {D, h, /3). (6.1) 

V ||r ^\\hs \\T ^\\hs I detT|/ 

Our task is to show that T is a scalar multiple of an orthogonal matrix. 

We may assume T — (*o 1^) is diagonal with ti,t2 > 0, by applying the 
singular value decomposition exactly as in the previous section. If ti = t2 
then the diagonal matrix T is a scalar multiple of the identity, and so the 
original T was a scalar multiple of an orthogonal matrix, as desired. 

Thus we may suppose from now on that ti ^ t2- We will deduce a 
contradiction, thus ruling out this possibility and completing the proof of 
the equality case in the theorem. 

Write u = ui for a first eigenfunction on D, with respect to the vector 
potential F{x) = ^{—X2,xi). Then u satisfies the eigenequation {ihV + 
F)^u = Xi{D, h, I3)u, which expands to say 

-h'^{uxixi+Ux2X2) + if>'l3{-X2,xi)-Vu + — {xj+xl)u = Xi{D,h,P)u. (6.2) 

Inspecting the proof of the theorem in the previous section, we see that one 
of the trial functions on T{D) is v = u o T^^, in other words v{xi,X2) = 
u{xi/ti,X2/t2)- (This trial function arises when m = N, since Un = Id.) 
Since equality must hold in the Rayleigh principle (j5.ip with n = 1, we 
deduce that this trial function must be a first eigenfunction on T{D) us- 
ing Planck constant bh and field parameter cf3. That is, v satisfies the 
eigenequation 

-b^h^{vxixr+Vx2X2)+i'bchl3{-X2,xi)-Vv^ — ■^{x\+x1)v = Xi{T{D),bh,cP)v. 



In terms of u, and using equality of the eigenvalues in (j6.ip . we can rewrite 
this last equation as 

2h / Uxixi , Ux2X2 \ I 2ih/3 . , \ ( Ux2 

+7;77JX72t(*?^? + *2^2)^ = Ai(D,n,/3)n. (6.3) 
i{ti + 12) 
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By solving the simultaneous linear equations (16. 2p and (16. 3p (which can be 
done since ti ^ t2) we obtain expressions for the second partial derivatives 
in terms of lower order derivatives: 

_ i/3 /32 2 , Xi{D,K/3) 

-Jl 2 ^ Ai(AM) 

Ux2X2 — ^XiUx2 4/j2^1^~'~ 2h^ ^' 

These formulas are equivalent to 

_^^iPxvx2/2h ^^^^^^ ^ ^2^^if)x^X2/2h ^g_4^ 
_^^iPx,X2/2h--^^^^^ ^^2(gi/3xiX2/2n-)^ (g^5) 

where w = ^/Xi{D,h, P)/2h^. Note w > because the eigenvalue Ai is 
positive (Lemma lA.Sp . 

Given any point y €z D, take a rectangular neighborhood R C D centered 
at y with the sides of the rectangle being parallel to the axes. On that 
neighborhood we can solve the ODEs (16. 4p and ()6.5p . finding that 

^il3x,x2/2h^ = a(x2)e^'^^i + 6(x2)e-^'^^S 
g-i/3xiX2/2rz,^ = c(xi)e^'^^2 + d(xl)e-^'^^^ 
for some functions a, b, c, d. Hence 

a(x2)e*'^^'i + 6(x2)e-*'^^i = e^'^^i^^Aj^^^^^^gic.xa ^ d{xi)e-'^''^) . (6.6) 

Let < zi < vr/cj with 4zi being smaller than the width of R. Con- 
sider equation ()6.6p with two different xi-values, namely xi = yi + z\ and 
x\ = y\ -\- 2zi, and with any X2-value that cuts through the rectangle R. 
The resulting two equations provide simultaneous linear equations for the 
unknowns a(x2) and 6(3:2). These simultaneous equations can be solved 
(since < 2ujzi < 27r). The solution for 0(^2) has the form 

0(^2) = e^l^y^^"^/^ X (hnear combination of terms e*'^^'^ with 

frequencies 7 = ita; + vjizi/h, where 1/ = 1, 2). (6.7) 

The same conclusion holds for any other sufficiently small value of zi , which 
implies that a{x2) is identically zero, as we now explain. 

Since /3 7^ we may choose d = ujh/\l3\k where the integer A; > 9 is 
taken so large that zi = 46 is smaller than tt/uj and 4zi is less than the 
width of the rectangle R. Let e = /3/\/3\ so that e = ±1. Then for the 
frequency 7 = ±uj + i^(3zi/hwe have e^"'^^ = e*(='='=+^^)('^^2/'=) when = I, and 
gt^x2 _ ^i{±k+8e){ujx2/k) ^j^gj-^ _ 2; these four exponentials determine the 
linear combination for 0(^2) in (j6.7p . On the other hand, we could choose 
zi = 5, in which case the linear combination for a{x2) would involve the 
four exponentials e^(±*^+=)('^^2/'=) and e^(±'=+2^)('^^2/'=). These two different 
linear combinations for a{x2) involve eight distinct frequencies: the eight 
integers ±k + e, ±k + 2e, ±k + 4e, ±k + 8e are distinct because k > 9. The 
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linear combinations for a(x2) are equal for a whole interval of X2-values, 
and so the coefficient of each of the eight exponentials must be 0, because 
a nontrivial trigonometric polynomial can have at most finitely many zeros. 
Hence a{x2) = 0. 

Similarly b{x2) must vanish identically and so u is identically zero on R, 
and hence on all of D, which is impossible because u is an eigenfunction. 
This contradiction concludes the proof. 

7. Proof of Corollary [m] 

Assume = T{D) for some linear transformation T. Pulling out the 
coefficient of the Planck constant in Theorem 13.11 implies that 

II IIH6 j=i I I j=i 

The Hilbert-Schmidt norm of T^^ can be interpreted in terms of moment 
of inertia of the image domain T{D), with 

Ik Whs ^ ' 
by [231 Lemma 5.3]. And obviously 

1 A{D) 

|detr| A{T{D))' 

Substituting these expressions now proves the desired inequality in the 
Corollary, after replacing (3 with I3/A{D) on both sides. 

The "maximality" statement follows from the "equality" statement in 
Theorem [O 

8. Proof of Theorem [33] 

The proof follows the Dirichlet and Neumann cases (see the proof of The- 
orem [3TTJ, except that we must handle also the boundary integral appearing 
in the Rayleigh quotient ([3.ip . This we now do, by following Laugesen and 
Siudeja's proof from the case of vanishing magnetic field [23|. Theorem 3.3]. 

The boundary contribution to the Rayleigh quotient of the trial function 
V on the domain E = T{D) is 

!QE\y?ds{x) ^ jQ^\u{UT-^x)\'ds{x) 
J^\v\'^dx J^\u{UT^'^x)\'^ dx 

_ jQ^\u{Ux)\''\Tr{x)\ds{x) 
J^\u{Ux)\'^dx-\deiT\ 

by X I—)- Tx, where t{x) denotes the unit tangent vector to dD at x. Geo- 
metrically, |rr(x)| is the factor by which T stretches the tangent direction 
to dD at X. 
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The rotational symmetry of D ensures that tangent vectors rotate accord- 
ing to t{U^^x) = U^^t{x). Thus after replacing x with U~^x in the last 
formula we obtain 

/g lyP^g(x) f \u(^^-j\2\TU-W{x)\ds{x). 

Jij 1^1 ox Jqd 

Choosing U = Um and averaging the preceding quantity over m, and then 
applying Cauchy-Schwarz, gives the upper estimate 

N 

{X] 



|detr|-i [ \n{x)\'{^ \TU-W{x)\'y' ds{, 

•'^^ m=l 

= |detr|-^ ^||r||H5 / \u{x)fds{x) 
v2 JdD 

by Lemma [4.31 After multiplying by aV^/\\T^^\\fjs aiid sustituting the 
relation 

\\T-'\\hs = j^^WTWhs 

(which is valid for 2x2 matrices, by the formula for T^^ in terms of T) we 
obtain the quantity 



fj / \u{x)\ ds{x). 
JdD 

Now we can easily modify the proof of Theorem 13.11 to handle the Rayleigh 
quotient ()3.ip for the Robin-de Gennes eigenvalues. 

The equality statement follows exactly as for Theorem 13.11 
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Appendix A. Invariance and positivity of the spectrum 

Our proofs rely on rotations, reflections and dilations of the domain. 
These reductions are justified by the invariance lemmas in this appendix, 
which we prove for the sake of readers unfamiliar with the magnetic Lapla- 
cian. 

Throughout this appendix, is a bounded plane domain, further as- 
sumed to have Lipschitz boundary in the Neumann and Robin-de Gennes 
situations. We give proofs for the Dirichlet eigenvalues only. The proofs 
are identical for Neumann conditions, except using instead of Hq. For 
Robin-de Gennes conditions, we indicate (when needed) how to treat the 
boundary integral in the Rayleigh quotient. 
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In the first lemma, we go against the notation used in the rest of the paper 
and write the eigenvalues as depending on the vector potential rather than 
on its curl, the mag netic field. Let F,F G C°°{R^;M.'^) be vector potentials 
defined on the whole plane that generate magnetic fields B = V x F, B = 
V X F, respectively. These fields are not necessarily constant. 

Lemma A.l (Gauge invariance). If two vector potentials differ by a gradient 
vector, then they generate the same magnetic eigenvalues. That is, if F = 
F + V/ for some function f G C°°{M?) then 

Xjin,h,F) = Xjin,h,F), j = 1,2,3,... 

and similarly for the Robin-de Gennes and Neumann eigenvalues. 

Proof. For the Dirichlet spectrum, given a trial function u G i?Q(i7;C) we 
define v = e^^ u G Hl{Q.;C). Then \v\'^ = |np and 

{ihV + F)v = e'f {ihV + F)u 

by direct calculation. Taking the magnitude and then squaring and inte- 
grating shows that the Rayleigh quotient of v with potential F equals the 
Rayleigh quotient of u with potential F. Hence the lemma follows from the 
variational characterization of the eigenvalues [3l p. 97]. □ 

Next we show that the eigenvalues depend only on the curl of the poten- 
tial, so that we are justified in writing \j{Vt, h, B) instead of Aj(il, h, F). 

Lemma A. 2 (Independence from choice of potential). If B = B on then 
F and F generate the same magnetic eigenvalues. 

Proof. Since is simply connected and the vector potentials F and F have 
the same curl, the potentials differ by a gradient vector. Hence they generate 
the same magnetic spectrum, by Lemma lA.ll □ 

Lemma A. 3 (Independence from direction of magnetic field). The fields B 
and —B generate the same magnetic eigenvalues. 

Proof. Given a trial function u to be used for the field B, define a trial 
function v = u ior the field —B. Obviously \u\ = \v\ and \{ihV + F)u\ = 
{{iKV — F)v\, and so the respective Rayleigh quotients of u and v are equal, 
which proves the lemma by the variational characterization of eigenvalues. 

□ 

For the remainder of the appendix we assume the magnetic field is con- 
stant, with B = (0, 0, /3) for some /3 e R. 

Lemma A. 4 (Invariance under rotation). If V is a 2 x 2 rotation matrix 
then Xj{V{Q,), h, /3) = Xj{Q,h,l3) for each j, and similarly for the Robin-de 
Gennes and Neumann eigenvalues. 
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Proof. Choose the vector potential F{x) = /3(Mx)^ = ^{—X2,xi) where 
M = ^ -^°2 ) 1 which gives the field V x F = (0, 0, /3). We use this vector 

potential on both the domains and V{Q). 

Given a trial function u on Q, define v = uo on Obviously u 

and V have the same norms, and thus the denominators of their respective 
Rayleigh quotients are equal. The numerator of the Rayleigh quotient for v 
has integrand 

\{ihV + F)v{x)\^ = \ih{Vu){V-^x)V-'^ + u{V''^x)F{x)\'^ 
= \ih{Vu){V-'^x) + u{V-^x)F{V-^x)\^ 

because F(x) = F(V~^x)V~^ (using that rotations commute in 2 dimen- 
sions and so commutes with V). Hence after integrating and changing 
variable with x i— )• Vx we see that the numerators of the Rayleigh quotients 
of V and u are equal. (Also, in the Robin case /gy(f^) ds = Jq^ |np ds.) 

The lemma follows immediately from the variational characterization of 
the eigenvalues. □ 

Lemma A. 5 (Invariance under reflection). IfVisa2x2 reflection matrix 
then Xj{V{i}),h, 13) = Xj{Q,h,(3) for each j, and similarly for the Robin-de 
Gennes and Neumann eigenvalues. 

Proof. In view of the rotation invariance in the preceding lemma, we may 
assume the reflection occurs across the X2-axis, with V = {~q i)- Choose 
vector potential F{x) = /3(0, xi), for which V x F = (0,0, /3), and check 
that F{x) = —F(V~^x)V~^ . Now adapt the proof of Lemma lA.41 to show 
Xj{V{il.),h, 13) = \j{Vt,h,—j3). Then replace — /3 with (3 on the right side, 
by Lemma lA.31 □ 

Lemma A. 6 (Invariance under translation). Ify G then \j{Q,+y, h, f3) = 
Aj(il, h, (3) for each j, and similarly for the Robin-de Gennes and Neumann 
eigenvalues. 

Proof. Straightforward, by translating the potential along with the domain. 

□ 

Lemma A. 7 (Invariance under dilation). If r > then Xj{ril.,rh, (3/r) = 
Xj{Q,h,f3) for each j, and similarly for the Neumann eigenvalues. For the 
Robin-de Gennes eigenvalues, pj{rQ,rh, l3/r,ra) = pj{n,h, I3,a). 

Proof. We use the potential F{x) = /3{0,xi) on il, giving held strength f3, 
and potential F{x/r) = 7(0, xi) on rf], giving held strength /3/r. 

Given a trial function u on Q., define v{x) = r~^u{x/r) on rO,. Then u 
and V have the same norms, and so the denominators of their Rayleigh 
quotients are the same. The numerator of the Rayleigh quotient for v is 



\{irhV + F{x/r))v{x)\''^ dx = I \ih{'Vu){x) + F{x)u{x)\'^ dx, 
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which equals the numerator of the Rayleigh quotient for u. In the Robin case, 
note also that the boundary term transforms according to ra Jq^^.^^^ \ v\'^ ds = 

Once again the lemma follows from the variational characterization of the 
eigenvalues. □ 

Finally we show the eigenvalues are positive, when the field is nonzero. 

Lemma A. 8 (Positivity of the energy). If f3 ^ then the first eigenvalues 
Xi{^l,h, /3), pi{il,h, f3) and fii{Cl, h, /3) are all positive. 

Proof. Suppose one of the first eigenvalues equals 0, and let n be a cor- 
responding eigenfunction using the potential F{x) = (0, Since the 
Rayleigh quotient of u equals 0, we deduce ihVu + Fu = in il. 

Consider a point at which u ^ 0. Near that point we have ikVlogu = 
—F, so that F is locally a gradient vector. Therefore F must satisfy the 
compatibility condition diF2 = d2Fi, which says /3 = 0. □ 

The first Dirichlet and Robin eigenvalues remain positive even when /3 = 
0. The Neumann eigenvalue /ii equals 0, when /3 = 0. 

Incidentally, monotonicity of the first Neumann eigenvalue ^i{il.,h,/3) 
with respect to the field strength is relevant in superconductivity. Mono- 
tonicity has been proved for large /3 by Fournais and Helff'er [9]. 
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